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Abstract 



O ■ It is known that the correlator of one axial and two vector currents, that receives leading contri- 

O: 



butions through one-loop fermion triangle diagrams, is not modified by QCD radiative corrections 
at two loops. It was suggested that this non-renormalization of the (VVA) correlator persists in 
' (— I '. higher orders in perturbative QCD as well. To check this assertion, we compute the three-loop 

QCD corrections to the (VAA) correlator using the technique of asymptotic expansions. We find 



^ _ that these corrections do not vanish and that they are proportional to the QCD /3-function. 



1. Introduction 

> 

■ The correlator of two vector currents and one axial current is an interesting object. Early 

OO ■ studies of this correlator were important for developing an understanding that, in spite of naive 

, equations of motion, the axial current is anomalous and that the perturbative part of the anomaly 

can be computed exactly from the one- loop triangle diagram The anomaly corresponds to 
just one term in the decomposition of the correlator into independent Lorentz structures, and the 
non-renormalization of the anomaly by higher-order QCD corrections does not say if other terms 
share the same property. In fact, since the discovery of the non-renormalization of the anomaly, 
it was generally believed that coefficients of other Lorentz structures change in higher orders of 
perturbation theory. This belief was challenged in Ref. |^, where it was pointed out that in 
the kinematic limit where momentum of one of the vector currents is vanishingly small, another 
non-renormalization theorem is valid. Indeed, in that limit just two independent form factors are 
needed to fully describe the (VVA) correlator. One of these form factors is the axial anomaly 
and, therefore, it is not renormalized. It was shown in Ref. 2] that due to helicity conservation 
in massless QCD, the two form factors are in fact proportional to each other, and so the non- 
renormalization of one of them implies the non-renormalization of the other. Therefore, we come 
to the conclusion that the entire (VVA) correlator is not renormalized in that limit. 

This result initiated new studies of the (VVA) correlator, which led to a number of surprising 
findings. First, in Ref. Q| it was pointed out that additional non-renormalization theorems for form 
factors of the (VVA) correlator exist even in the case of the most general kinematics. Later, in 



Ref. 



an explicit calculation of the two-loop 0{as) contribution to the correlator was reported. 
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The result turned out to be extremely simple: it was found that the 0{as) correction to the 
correlator vanishes identically in the most general kinematics if a consistent definition of the axial 
current is employed. This unusual feature of radiative corrections prompted the authors of Ref. 
to speculate that their result might be an early indication of the non-renormalization of the entire 
(VVA) correlator in perturbative QCD with massless quarks to all orders in the strong coupling 
constant. 

The goal of this Letter is to show that this speculation is not correct and that radiative 
corrections to the (VVA) correlator appear at the three-loop order in perturbative QCD. These 
three-loop corrections are, however, peculiar in that they are explicitly proportional to the QCD /3- 
function and vanish in the conformal limit of QCD. This is a natural result. In fact, it was pointed 
long ago in Ref. [d] that the exact (VVA) correlator in a conformally-invariant theory is given by 
the one-loop expression and no contributions at higher-loops are allowed. This happens because 
conformal symmetry restricts the functional form of the three-point function up to a possible 
multiplicative factor which, however, is fixed by the requirement that the non-renormalization of 
the anomaly works out correctly. It is easy to realize that through two loops in perturbative QCD, 
no diagrams that cause violation of the conformal symmetry contribute to (VVA) correlator, while 
at three loops diagrams that correspond to the running of the coupling constant explicitly appear. 
This allows us to understand why no corrections to (VVA) correlator were found in Ref. ^, 
why such corrections appear at the three-loop order, and why they are proportional to the QCD 
/3-function. 

The remainder of this Letter is organized as follows. In the next Section we define the (VVA) 
correlator and discuss its decomposition in terms of invariant form factors. In Section[3]we explain 
how three-loop contributions to the correlator can be computed. In Section|3]we present the results 
of the calculation. We conclude in Section [5] 

2. Definitions 

We consider the correlator of two vector currents and one axial current 

W^,,p{qi,q2) = - /'dVd4x2e'(9--i+9--^) ( | T{l^^(a;i)K(x2)Ap(0)} | ). (1) 



2 



The currents are flavor-diagonalQ They read 

Vf, ^ ipif^i} , ^ ipltils ip , (2) 

where "0 are fermion fields charged under the SU{Nc) color group. 

The {VVA) correlator in Eq.([T]) can be expressed in terms of independent Lorentz structures 
and invariant form factors Q . A suitable parameterization reads 



W^.p(gi,<?2) = WL{qlqlq'')ti'^p + W^+\qlqlq%+Jp 

+ w'f^ {ql , ql , W + w^f\ql , qj , , 

where, following Ref. we introduced four independent tensor structures 

'■/xj/p ~ qp^fj,iyal3qi </2 I 

tjZp = qiu^tLpami I2 - q2^,.^upapqi q2 -qi- q2^^,.ypaK - 2 e^t^affgi 92 9P ■ 

f(-) _ a 13 _ ^ ' 1 a P 

pvp — £/^i/Q/3'?l 92 '^P ^2 ^^"^"(991 92 9p ' 



(3) 



(4) 



tliup = qiuefj.pai3qiq2 + q2t,ei'pai3q?q2 - qi ■ q2etj.i^paq°' , 

and momenta q = qi+q2, n = qi—q2- From momentum conservation, q is the incoming momentum 
carried by the axial current. We note that all tensor structures are transversal with respect to the 
momenta of the vector currents 

qit^^p = 0, q2tp„p = , (5) 

where t^^p is a generic notation for t'jjjp, t^^\ and ijitp. On the other hand, these tensor structures 
have different transversality properties with respect to the momentum q of the axial current, 

q''t^pip = ~q^e,,o.pqUL q't^^l^^, q'Vp ^ ^ ■ (6) 

Therefore, the parameterization in Eq.Q is consistent with the conservation of the vector current 
and with the fact that the conservation of the axial current is violated by triangle loop diagrams. 
We will refer to wi, as the longitudinal form factor and to w!^"' and w!^ ' as the transversal form 
factors. 

We now discuss what is known about these form factors. In general, they are functions of 
three independent kinematic variables, gj, g| and q^ = {qi +92)^- Since the divergence of the axial 
current is fixed by the anomaly equation 

d.A'^-^e.^.^F'^^FP'^, (7) 



^Note, however, that wc do not include the so-called singlet contributions to the correlator when we compute it 
in perturbation theory. 
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to all orders in perturbation theory, the longitudinal form factor is completely defined. Indeed, 
the anomaly equation implies 



27r2 



We solve this equation for wl to find 



/ 2 2 2 

WL{qi,q2,q 



iNr 



(8) 



(9) 



27r2g2 • 

While the longitudinal form factor can be fully determined from the anomaly equation, the 
transversal form factors can not. However, as it was shown in Refs. , some combinations 
of transversal form factors are uniquely fixed by the chiral symmetry in perturbative QCD with 
massless fiavors. The following relations are valid jsj. 



,„(+) 



■ w, 



(-) 



.,(-) 



r,(-) 



- w. 



2 2 2\ 

iql,ql,q^) + 
{qlqlq^) + 



.,(+) 



- W, 



{q\ql,ql)^Q, 
{q\ql,ql)=Q, 

(q ,92>'7i) - WL[q ,g2,9i) 



(10) 



,92 •9„„( + ) 



qf 



);>{q\qi,qi)-2 



qi ■ '72 (-) 



qf 



w'r>{q',qlqt) 



In perturbation theory, the correlator W^iyp receives contributions starting at one-loop order, 
see Fig. [TJ The correlator was computed through two loops in Ref. 0] for arbitrary momenta 
qi and q2. The results of that calculation showed that, with the proper definition of the axial 
current, the two-loop corrections to the full correlator vanish. We pointed out a possible reason 
for that in Section 1, where we noted that it is natural to expect that these results will be violated 
starting from three loops. In the next Section we discuss how to set up a calculation to check this 
explicitly. 



3. The three-loop calculation 

We are interested in computing the Wp^p correlator to three loops in perturbative QCD, for 
arbitrary qi and q2. Sample diagrams are shown in Fig.[TJ We note that the one- loop contribution is 
described by 2 diagrams, the two-loop contribution by 12 diagrams and the three-loop contribution 
by 182 diagrams. It goes without saying that, at present, no computational technology exists that 
can be used to calculate three-loop three-point functions with three external off-shell legs. To 
circumvent this problem, we consider a hierarchy of incoming momenta, 52 ^ '7i, and construct 
a systematic expansion of the three- loop diagrams in the ratio 92 /^i- It is intuitively clear that, 
upon such an expansion, a three-point function is mapped onto a set of two-point functions and 
its derivatives evaluated either for small k ^ q2 or large k ^ qi loop momenta. Such an expansion 
is constructed most efficiently by using standard techniques of the large momentum expansion ^ . 
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In the current application, we use this procedure as automated in the computer programs q2e 
and exp [lo| . Once the three-loop three-point functions are mapped onto the three-loo p tw o- 
point functions by the large momentum expansion procedure, we use the program MINCER [12] to 
compute them. Finally, we note that all diagrams that contribute to the correlator are generated 
with Qgraf Q. 

As a technical remark, we note that treatment of tensor integrals that appear in this computa- 
tion requires some care. At one loop, MINCER routines can deal with tensor integrals of arbitrary 
rank, contracted with arbitrary external momenta. Unfortunately, at higher loops MINCER can 
only be used to compute two-point scalar integrals that depend on a single external momentum. 
However, since we arrive at relevant two-point functions by expanding three-point functions in 
q2/qi, we need to compute more complicated tensor integrals. To give an example, consider an 
integral 

^ d'^fc, (fci-g2)°^(fc2-g2)°^. , . 

I\i27rrklkl..{ki-qir..{k2-qir..' ^ ' 

where the two-point function itself depends on momentum qi but the numerator contains scalar 

products of loop momenta ki with another external momentum 52 • This is a typical situation for 

us to face, and it can not be handled by MINCER. To deal with the integrals of the type shown in 

Eq. ()lip . we write qi^ — (92 • 9i)/'?i'?f + 92 _L' where 52, _l • 9i = 0, and note that the final result can 

not depend on the direction of 92. We make use of this observation by averaging over directions 

of q2.i_ in Eq. (jll[) . Such averages are constructed from the metric tensor g'^" — q^q^ /q\ which only 

contains vector qi and, therefore, leads to scalar products of four-momenta that can be treated by 

MINCER. 

Computational procedures described above rely heavily on the use of dimensional regulariza- 
tion. This leads to a subtlety since, as it is well known, continuation of the Dirac matrix 75 to d 
dimensions requires care and a proper way to do so is crucial for the correct computation of the 
lyVA) correlator. We use the definition of the axial current developed in Ref. Q which, on one 
hand, is self-consistent and, on the other hand, does not introduce unnecessary difficulties into 
already complicated multi-loop computation. Thus, we define the axial current in the following 
way [el, 

Ap = '^ep^,p,x4T^^'P'^'H, with £0123 = 1 = -e"'"' , (12) 
where T"^^ = j'^'-fl^j'^ and square brackets indicate the anti-symmetrization of the Lorentz indices, 
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With this definition, the axial current needs renormalization. We write Ap — Za^^^° and note 
that the renormalization constant Za was worked out in Ref. i6;]. There, in addition to the MS 
renormalization constant, a finite piece is added that ensures that the non-singlet current has no 
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anomalous dimension and the anomaly equation ([7]) is satisfied. The full renormalization constant 
reads at two loops 

Za = (^1 + ^^^/3o) x(l- 4CFa, + Cpaj (^22 Cp - ^ Ca + ^Uf^^ , (14) 

where = g^(/i)/(167r^) is the MS QCD coupling constant defined at the scale fi and /3o = 
11/3Ca — 2/inf is the one-loop QCD /?- function. Also, n/ is the number of (massless) quark 
flavors and Cp — {N^ — l)/(2iVc), Ca — Nc are the Casimir operators of the SU(Af) group in 
fundamental and adjoint representations, respectively. In QCD, we use iV = 3. 

We now return to the definition of the axial current in Eq. ([T^ and observe that it contains the 
Levi-Civita tensor e^^ap, which is a four-dimensional object whose continuation to d dimensions is 
not possible. To address this issue, we first compute an auxiliary tensor correlator Tai3\,p.v which 
does not contain the Lcvi-Civita tensor and which is finite at d = 4. We then obtain the required 
correlator by contracting Tai}\,^v with the Levi-Civita tensor 

y^p.vp = -^^paPxTapX^^iv ■ (15) 

This can be easily done because all entries on the right hand side of Eg. p^ are well-defined in 
four dimensions. 

To define the auxiliary tensor Ta^x.^u, we imagine that the correlator W^^p is contracted with 
^apxp 

Tai3x,f,i^{qi,q2) = -eQ^ApVVA'i/p(<7i,g2) ■ (16) 

We then use the identity 

^a/lXpipa'P'X' = ga,a'gi3,l3'gX,X' — Qa^a' 9 ,X' gX,P' + • ■ • (17) 

to remove the Levi-Civita tensors in favor of products of metric tensors and find a convenient 
expression for the auxiliary correlator 

Tal3X,p.u{qi,q2) {9aa'90l}'9XX' " 9aa' 9 13X' 9XP' + ■■) 

r , , , (18) 

X / d^xid'^X2e'^'''-'''+'i^-''^^ (0|T{l/p(xi)K(a;2)^(0)r["'^^lV(0)}|0). 



Because the right-hand side of Eq. (|18p can be analytically continued to d dimensions in a straight- 
forward way, the correlator Tapx.fj.uiqi, q2) can be calculated in dimensional regularization. After 
the renormalization of the "axial" current Ea. (jl4l) is applied, finite result for Tai3x.ij.v{qi,q2) is 
obtained. Finally, we use Eg. dTSl) to calculate the correlator W^i/p. 

Before discussing the results of the computation, it is useful to comment on the Lorentz de- 
composition of the correlator Tapx,p,u- As follows from its definition, Eq. ()18p . Tai3x,tiv{qi,q2) is 
an anti-symmetric tensor with respect to its first three indices and it is a symmetric tensor under 
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a simultaneous change o 52, M ^ ^- We find that it is possible to express T in terms of six 
tensor structures 

6 

Tapx,^,u{qi,q2) = X! "=^(91' 92) T''alx,t,v ' (19) 
1=1 

where 

The square brackets here imply an anti-symmetrization with respect to indices a, /3 and A while 
the positions of /i and i' are kept fixed. 

While the tensor Tai3x,fii^ can be conveniently computed using the asymptotic expansion tech- 
nique, some care may be required with the interpretation of the result, since the coefficients 
ci-el'Zij 92) are not independent. Relations between them can be obtained from the conservation 
of the vector currents and from the non-renormalization of the axial anomaly. We find three 
relations 

C2 - 91 C4 - (71 • 92 C6 = , ci - gi • 92 C3 - C5 = , a - C2 = ^ ■ (21) 

In the actual three- loop computation, Eqs. ipTj) are not enforced but are used as checks of its 
correctness. Once coefficients ci„6 are known, we compute Wapx using Eq. ljlSp and re-write it 
through the form factors if^^^ ■* and wl- The corresponding results are presented in the next 
Section. 

4. Results 

We are now in position to present the results of the calculation. To this end, we introduce the 
notation 

qi-q-2. I2 T 1 ^oo^ 

ri = —, r2 = ^, Lr^lnr2. (22) 

qi qf 

We perform the calculation in the limit g2/<Zi ^ 1, so that ~ r2 ^ 1. The results for the form 
factors below are given as an expansion in ri_2- For each transversal form factor we find that 

• the two- loop QCD corrections vanish, in accord with the observation of Ref. ^j; 

• the three-loop QCD corrections do not vanish, but arc proportional to the one-loop QCD 
/3-function0 As a result, they vanish in the conformal /3 — ;> limit. 



^We note that the simplest way to check that the three-loop corrections do not vanish is to study nj--dependent 
vacuum polarization insertion diagrams. Since the correction to the correlator turns out to be proportional to the 
/9-function, such diagrams give, in fact, the full answer. 
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The longitudinal form factor is known exactly from the anomaly equation; it is given in Eq. 
We do not discuss it anymore. The perturbative expansion for a generic transversal form factor 
wt — , is written as 

WTiqlqlq') = (wi^:^ + PoCf (^^J J^^) + 0{a^)^ , (23) 

f 1 3) 

where /3o — 11/3Ca ~ 2/3n^. We calculated wj^' using asymptotic expansions in (72/91 through 
fourth order. Upon expanding the one-loop expressions for these form factors presented in Ref. [J|, 
we find full agreement with our result |f| 

Below we present the expansions of the one-loop and three-loop form factors through second 
order in q2/qi- For the one- loop contributions we find 

(+ 1) 20 4 ^ 2ri , ^ , 2/^52 8 \ / 82 8 \ ^^,3, 

(-!■) 16 4 2ri , , ^ /84 8 \ /48 4 \ , /../i^ 

= -«,(r'^' + 0{rl) . 

The three-loop contributions read 

4+^^) = 2L2 + (-14 + 16C3)ir + y - f C3 - ri (^2^2 - (10 - 16C3)i. + (^^ " ^ 

2 /3853 128^ /91 96, \ ^ 22 . 
+ '^^ + — C3 - I — - — C3 ) + — 



V 270 25 ^ V 9 ^ J 9 
/ 6991 88 , /125 32 \ 11 ,\ 

4-'^) = 2Ll + (-14 + 16C3)i. + f - f C3 - (2L^ + (-10 + 16(3)^. + f - ^(3 ) (^5) 

2 /4573 128 /91 96 \ 22 . 



f? + Cs Cs Lr + —L: 

^ V 270 25 V 9 5 ^7 9 

/^9209 496, /^127 48, \ ^ ,7^2 



-(-.3) (-3) 

= — 



0(r?) 



It is now straightforward to check the exact perturbative QCD relations between different form 
factors, shown in Eg. dTU]) . The above expansions are given for wriqi, q2, q^)] to compute Eq.([TU]), 
we require wriq'^, 92i These can be easily obtained from the above equations by applying the 
following transformations to them 

ql^ql{l + 2n+r2), ^ r^/ {I + 2n + r2), n^- '^^^ . (26) 

1 -I- 2ri -I- 7-2 



^We note that this agreement is only found if we multiply all results in Ref. 3 by a factor two that, we believe, 
is erroneously omitted there. 
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Upon applying these transformations to Eas. (j24l25l) . re-expanding in ^ r2 <^ 1 and combining 
the form factors in the right way, we find that Eqs.([TU|) are indeed satisfied. As a simple illustration 
of this fact, we note 3] that, in the limit (72 — ^ 0, these equations reduce to 

WLiqlO,ql)^2{w^+\ql,0,ql)+W^f\qlO,ql)) . (27) 

Since the longitudinal form factor does not receive QCD corrections, the three-loop corrections to 
w^-* and w!^ ^ in these kinematics must cancel. Taking the limit r2 — > 0,ri — > in Ea. (|25p we 
observe the required cancellation. 

5. Conclusions 

In this Letter, we described the calculation of the three-loop perturbative QCD contribution 
to the correlator of one axial and two vector currents for general kinematics. In perturbative QCD 
this correlator receives contributions starting at the one-loop order. The two-loop contribution 
was computed in Ref. 4|, where it was observed that, for a propely defined axial current, the 
two-loop contribution vanishes. This result prompted the authors of Ref. to suggest that this 
non-renormalization of the full correlator will carry over to yet higher orders in perturbative QCD. 

Our three-loop computation clarifies this issue. We compute the three-loop contributions to the 
correlator in arbitrary kinematics, using the techniques of asymptotic expansions. We find that the 
three-loop contributions do not vanish but that they are proportional to the QCD /3-function. This 
feature explains the vanishing of the two-loop contribution to (VVA) as being due to conformal 
symmetry of QCD. Indeed, it was pointed out in Ref. [5] that in the conformally-invariant theory 
the functional form of the correlator is fixed. Through two-loops the perturbative contributions 
to (VVA) in perturbative QCD are not sensitive to the breaking of conformal symmetry, but 
this changes at three loops because diagrams that describe the running of the coupling constant 
appear for the first time. Thus, the three-loop contribution should be proportional to the /?- 
function, in full accord with our explicit computation. Finally, we have verified the validity of all 
non-renormalization theorems in perturbative QCD for the form factors of the (VVA) correlator 
derived in Refs. 0, We find that through three- loops, all these theorems are satisfied. 
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